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ABSTRACT 

We examine the temperature structure of the intergalactic medium (IGM) surounding 
a hard radiation source, such as a Quasi-Stellar Object (QSO), as it responds to 
the onset of helium reionization by the source. We model the reionization using a 
radiative transfer (RT) code coupled to a particle-mesh (PM) A^-body code. Neutral 
hydrogen and helium are initially ionized by a starburst spectrum, which is allowed 
to gradually evolve into a power law spectrum (oc z/""'^). Multiple simulations were 
performed with different times for the onset and dominance of the hard spectrum, 
with onset redshifts ranging from z — 3.5 to 5.5. The source is placed in a high- 
density region to mimic the expected local environment of a QSO. Simulations with 
the source placed in a low-density environment were also performed as control cases to 
explore the role of the environment on the properties of the surrounding IGM. We find 
in both cases that the IGM temperature within the Hem region produced exceeds the 
IGM temperature before full helium reionization, resulting in a "thermal proximity 
effect," but that the temperature in the Hem region increases systematically with 
distance from the source. With time the temperature relaxes with a reduced spread as 
a function of impact parameter along neighbouring lines of sight, although the trend 
continues to persist until z = 2. Such a trend could be detected using the widths of 
intervening metal absorption systems using high resolution, high signal-to-noise ratio 
spectra. By contrast, the Doppler widths of Hi absorption features in mock spectra 
along neighbouring lines of sight show a weak trend with impact parameter prior to 
full helium reionization reflecting the behaviour of the underlying density and peculiar 
velocity fields, but they take on a near constant value after helium reionization, with 
a median value near 25 km s^^ at z = 3, in good agreement with observations. 

Key words: radiative transfer - quasars: absorption lines - quasars: general - cos- 
mology: large-scale structure of Universe - methods: N-body simulations 



1 INTRODUCTION 

An abundance of observational evidence shows that sub- 
sequent to the Recombination Era, during which both hy- 
drogen and helium recombined to their neutral forms, the 
Universe was reionized. The best constraints on the epoch 
of hydrogen reionization are derived from measurements of 
the intergalactic Lyman resonance line transitions in the 
spectra of high redshift Quasi-Stellar Objects (QSOs) and 
polarization measurements of the Cosmic Microwave Back- 
ground (CMB). The measured optical depths of intergalac- 
tic Lya and Ly/? absorption in high redsh ift QSOs show 
that the Universe wa s reionized by z ^ 6 (iBecker fc et all 
I2OOII : iFan et al.ll2002l '). This is consistent with CMB mea- 
surements by the Wilkinson Microwave Anisotropy Probe 
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(WMAP), which place the hydrogen reionization epoch, if 
a sudden process, at z = 11.0 ± 1.4 l|Dunklev et al J 120091 '). 
The sources of hydrogen reionization are currently unknown, 
but are strongly suspected of being young star-forming 
galaxies. More speculative possibilities include pockets of 
Popu lation III stars, as may arise in isolated stax clus- 
ters dMadau et al.l Il999l: [Choudhurv fc Ferraral l2006l h or 
miniquasars i Madau et al.ll2004 ). Hydrogen reionization by 
QSOs is untenable based on current estimates of the QSO 
luminosity function, unless there is a sharp rise at the faint 
end at high redshifts ( Meiksin 20051 : ICowie et al, 20091. 

The reionization epoch of helium, by contrast, is far 
less well constrained. This is primarily because of the need 
to go to space to measure the intergalactic Lyman series ab- 
sorption of intergalactic gas in the spectra of high redshift 
QSOs. The current constraints place the helium reionization 
epoch, when helium becom es nearly fully doubly ionized, at 
z >3 (jReimers et al.ll2005l ). Because of the high photoelec- 
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tic threshold energy of Hen, it is expected that helium was 
reionized by QSOs rather than stars, most likely in the red- 
shift interval z = 3 — 4 |Madau fc MeiksinI IT993 : iMeiksinI 
I2OO5I : iMcQuinn et alJlioog T 

Measurements of the Hen Lya optical depth show 



the optical depth rises rapidly appro aching 



3 with 



an increasing amo unt of patchiness (jZheng et al.l l2004l : 
iFechner et al. l2006f). By z ~ 3.4, t he optical depth is im- 
measurably large ( Zheng et al.ll200§ ) . The limits are not yet 
strong enough to demonstrate the epoch of helium has been 
detected, but the rapid rise in optical depth and degree of 
patchiness are suggestive of an approach to the helium reion- 
ization epoch. This interpretation, however, is not unique; it 
may also indicate a diminishing attenuation length of Heii- 
ionising photons, resulting in large fluctuations in the Hen- 
ionising metagalactic background due to large local Poisson 
fluctuations in the numbe rs of QSOs contributing to the 
background (|Meiksinll2067l ). The uncertainty in QSO counts 
and the spectral shape of high redshift QSOs still does not 
allow much ear lier reionization, up to 2 ;^ 5.5, to be excluded 
l|Meiksinll2005h . 

In this paper, we examine an alternative means of prob- 
ing the epoch when helium in the intergalactic medium 
(IGM) surrounding a QSO was reionized: detections of a 
heightened IGM temperature above the post-reionization 
equilibrium value. The temperature following full helium 
reionization will reach high values as the helium ioniza- 
tion front (I-front) passes, eventually establishing a lower 
temperature in ov erdense regions as the gas achieves ther- 
mal equihbrium (iMeiksinI 1 1994 iTittlev fc MeiksinI l2007l : 
iMcQuinn et~aLl |2009| ). If the IGM near the QSO is pol- 
luted with metal lines, then a direct determination of the 
gas temperature may be made by combining the measured 
linewid ths of the m etals, provided they are adequately re- 
solved (|Rauch et al . 1996; Meiksin 200"^). Because of the low 
abundances of intergalactic metals, however, measurements 
have so far been restricted to high Hi column density sys- 
tems. The advent of large telescopes, like the Thirty Meter 
Telescope^ or an Extremely Large Telescope!!!, may permit 
temperature measurements to be extended to much smaller 
column densities. 

Whilst such direct temperature estimates are currently 
limited to high column density absorbers, the reionization 
of helium by a QSO will induce a "thermal proximity ef- 
fect" in which the increased temperature near the QSO 
will be reflected by broadened hydrogen, as well as helium 
and metallic, absorption lines. The effect has been exploited 
as a means of probing the epoch of full helium reioniza- 
tion through the global statistical effect on the distribu- 
tion of the Hi Doppler parameters for systems throughout 
the IGM. Rises in the lower cutoff of the Doppler parame - 
ter distribution h ave been reported bv lRicotti et al.l (|2000l ) 
and lSchave et all ( 2000), and interpreted, based on calibra- 
tions with numerical simulations, as evidence for a sharp 
temperature jump and reduction in the effective polytropic 
index towards isothermality near z — 3.0, due to the on- 
set of_JIei^j;eioiiization. A similar conclusion was reached 
bv lZaldarriag^ ^M) 

in a wavelet-based analysis of QSO 
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spectra. Global reionization simulations and semi-analytic 
estimates allowing for distributed sources, however, are not 
clearly in accord with such sharp changes in the absorp- 
tion line parameters. The results suggest only a gradual 
rise in the mean I GM temperature should occur during he - 
lium reionization l|McQuinn et al.ll2009l : iBolton et al.ll2CI09l '). 
Searches for isolated regions of increased temperature, as 
may arise around a helium-ioni sing QSO, have been only 
partially successful. IZaldarriagal |2002l argues at least ten 
QSO sitelines are required to probe temperature fluctuations 
at a statistica lly signiflcant level. Using a similar method, 
iTheuns et al. flnd a region 10* km s"^ across sug- 

gestive of a statistically significant elevated temperature in 
a single QSO siteline out of eight examined. 

A more direct means of searching for the reionization 
signature is to target the region around specific QSOs. The 
temperature of the IGM gas surrounding a QSO may be 
probed both using lines of sight towards the QSO and nearby 
on the sky towards background QSOs. Surveys of QSO 
neighbours with follow-up spectroscopy along multiple lines 
of sight are currently underway in an attempt to measure the 
transverse hydrogen ionization proximity effect and to use 
transverse spatial correlations in the hydrogen Lya forest as 
a constraint on the vacuum energy density of the Universe 
(|Marble et al.l l2008t ) . The detection of a thermal proxim- 
ity effect would provide additional support for unification 
models of Active Galactic Nuclei (AGN). In these models, 
the variety of AGN activity observed, from radio galaxies to 
QSOs and from Seyfert Is to Seyfert 2s, is a consequence of 
geometry: an obscuring dusty torus or other obstacle blocks 
the optical/UV radiation from the central engine along some 
lines of sight, possibly re-radiating it in the infra-red. The 
observed character of the object then depends on the view- 
ing an gle. While good evidence exists for elements of the 
model ijVerma et al.ll2005l : lAntonuccil!l993l '). direct evidence 
for optical/UV radiation transverse to the line-of-sight of 
a putative obscured QSO or Seyfert 1 is generally lacking. 
Prior to the completion of helium reionization in the IGM, 
the broadening of Lya absorption features transverse to the 
line of sight toward bright radio or infra-red galaxies. Type 
II QSOs, or Seyfert 2 galaxies would provide evidence for in- 
tense photoionising UV radiation beamed transverse to the 
line of sight, as predicted by the unification models. 

The purpose of this paper is to use numerical reioniza- 
tion simulations to quantify the impact of Hen reionization 
by a QSO on the temperature and absorption line widths 
produced in the surrounding gas. The simulations are per- 
formed using a radiative transfer code coupled to an A^-body 
code to model the evol ution of the IGM. The m ethod was 
previously described in iTittlev fc MeiksinI (|2007l ). The de- 
tails of the simulation volume and sources are given in 
Results of the simulations are provided in fJ^IS'iid discussed in 
^ Our conclusions are summarised in ^ In an Appendix, 
we describe modifications to the radiative tranfer code im- 
plemented to simulate reionization by a central source. 



2 SIMULATIONS 

We use PMRT l|Tittlev fc MeiksinI 120071 ') to simulate the ex- 
panding ionization front produced by a central source. The 
simulation code is the merger of a Lagrangian particle-mesh 
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Figure 1. The evolution of the gas temperature surrounding a central radiation source which evolves from a starburst spectrum to QSO 
(power-law) over the redshift interval 4.5 > z > 3.5. The source is at a density maximum in the simulation volume. The images show the 
rapid increase in temperature following the turn on of the QSO. The polar grid radius is 25/2^/^ Mpc (comoving). 



(PM) code (|Meiksin et al.|[l999l ). and a grid-based radia- 
tive transfer (RT) co de based on a probabilistic photon - 
conserving algorithm (|Abel et al.lll999l : iBolton et al.ll2004 ). 
The method assumes the baryons trace the dark matter. 
The code is thus not able to cope with sudden increases 
in the pressure forces as the gas is heated and driven out 
of shallow potential wells or heated by sho cks in collapsing 
structures jBond et al.lll988l : lMeiksinHl994l '). Whilst pseudo- 
hydrodynamical methods are no replacement for full hydro- 
dynamics computation for high-precision comparisons with 
data, comparisons with full hydrodynamics computations 
show that the method recovers the distribution of Doppler 
parameters to high accuracy, although the median Doppler 
arameter may be slig htly broadened by 1 — 2 km s~^ 



Meiksin fc Whitdl200ll ). which is an intrinsic limitation to 
our estimates. An alternative method including an enthalpy- 
based pseudo-pressur e term similar to the HPM scheme of 
iGnedin fc HuJ (|l998l ) was found to broaden the lines even 
more. For this reason we adhere to the more s traightfor- 
ward implementation of lMeiksin fc White! l|200ll ). In future 
full hydrodynamics computations including radiative trans- 
fer would be preferable, and the authors are developing such 
a scheme. The pseudo-hydrodynamics scheme none the less 
does have advantages over a fully hydrodynamical one: it 
runs much more quickly without the hydrodynamics over- 
head, and the reduced memory requirements permit simula- 
tions in larger boxes at the high spatial resolution required 
to resolve the absorption systems. The latter is important to 
ensure the simulations both capture the large scale modes 
required to converge on the Doppler parameters and to re- 



solve the internal structure of the systems jTheuns et al.l 
1 19981 : iBrvan et al.lll999l : iMeiksin fc Whitell200ll ). 

The central QSO is modelled as a starburst that gradu- 
ally develops into a hard QSO spe ctrum. The star burst spec- 
trum was produced by PEGAS^ ( Fioc fc Rocca-Volrn orang^ 
I1997I ) for a galaxy 30 Myr after a burst of Pop III (zero 
metallicity) star formation. The starburst spectrum has a 
luminosity at the hydrogen Lyman edge of 10^* WHz^^. It 
is able to ionize hydrogen and some neutral helium to singly 
ionized helium, but has negligible flux above the ionization 
threshold of singly ionized helium. The spectral index of 
QSO sources at these high energies, particularly above the 
Hen Lyman edge which is most relevant for this work, is 
virtually unconstrained by observations. Using HST QSO 
spectra cove r ing p rimarily the redshift range 0.3 < z < 2.3, 
iTelfer et al.l l|2002l ') report a steep source spectral index of 
as ~ 1.76±0.12 for a fit over the spectral range 0.8—1.8 Ry, 
still a factor of two short of the Hen Lyman ed ge. By con- 
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trast using FUSE data for QSOs at z < 0.67, IScott et all 
(j 2004| ) find a spectral index of as = O.SOjgg^® fit over the 
similar spectral range 0.8 — 1.4 Ry. A fair fraction of the 
QSOs have quite hard spectra with as < 0. The origin of 
the discrepancy with the HST sample is unclear. The FUSE 
sample is somewhat fainter, although there is substantial 
overlap in the QSO luminosities with the HST sample. A 
correction for any intervening Lyman limit systems and a 
statistical correction for intervening intergalactic absorption 
are applied to the HST spectra, which is unnecessary for the 
lower redshift FUSE sample, so that the FUSE sample may 
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Table 1. Identifiers for tlie various source models used in this 
paper. 



Figure 2. Tlie density contrast 5p = {p — {p))/{p) at z = 2.6 
surrounding a central source positioned at a density maximum in 
the simulation volume. The shading ranges from 5p = —1 (black) 
to Sp > 100 (white). The polar grid radius is 25/2^/2 h'^ Mpc 
(comoving). 



be a truer reflection of the intrinsic QSO spectra at these 
energies. Ahernati vely, it may be that QSOs harden at lower 
redshift, although Pfelfer et all (|20o5) find no such effect in 
their data. Even if most QSO spectra tend to be softer at 
higher redshifts, for a given luminosity at the Hen Lyman 
edge QSOs with harder spectra will provide a greater num- 
ber of ionizing photons, and so they may still dominate the 
reionisation of helium in the IGM. The QSO spectrum in our 
simulations is approximated as an as = 0.5 power law at en- 
ergies above the Hi Lyman edge. Although a hard spectrum 
is suggested by the FUSE data, our choice is primarily mo- 
tivated by an attempt to address the discrepancy between 
the measured Hi Doppler parameters and the predictions 
of simulations. Detailed comparisons between simulations 
and data suggest the simul ations substantially underpredict 
the t emperature of the ga s (jTheuns et al.ll 19981 : ISrvan et al.l 
1 1999i : Meiksin et al.ii200l] ). A harder spectrum will produce 
a higher temperature. If still inadequate, then additional 
sources of heating or line-broadening (such as winds) may 
be required. Ultimately an assessment of the impact of Hen 
reionization on the gas temperatures will require improved 
knowledge of the intrinsic spectra of QSOs, particularly at 
energies above the Hen Lyman edge. 

The QSO spectrum is normalised to provide a peak 
hydrogen-ionising photon production rate half that of the 
starburst, with a peak luminosity at the hydrogen Lyman 
edge of 2.6 x 10^^ WHz~^. Three transitions from starburst 
spectrum to QSO were simulated: a gradual transition over 
the redshift intervals z = 5.5 to z = 4.5, z = 4.5 to z = 3.5 
and 2; = 3.5toz = 2.5. The starburst spectrum was turned 
on at z = 8 in all cases. In principle, an ionization front 
representing the accumulated ionising photons from a pop- 
ulation of galaxies could have been swept across the simu- 
lation volume, however we find that the IGM temperatures 
have relaxed by the time the QSO spectrum turns on, with- 
out any memory of the central source, so that initiating hy- 
drogen ionization by a central source adequately produces 
the expected temperature structure of the IGM prior to the 
turn-on redshift of the QSO spectrum. In post-processing 
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the spectra, we have added a uniform hydrogen ionization 
rate to match recent estimates. 

To approximate the expected environment of a QSO, we 
computed mass overdensities in the simulation after smooth- 
ing the PM distribution at z = 3.2 over a comoving scale 
of 0.5/i~^ Mpc, and placed the source at an overdensity 
peak corresponding to a mass of 1.5 x 10^^ Mq, consistent 
with recent estimates of QSO dark matter halo m asses based 
on QSO space density a nd clustering statistics dShen et al.l 
I2OO7I : IWhite et al.|[2008l ). In order to gauge the importance 
of the environment to the resulting properties of the IGM, 
we also computed a set of models with the source placed at 
the centre of an underdense region. Qualitative effects found 
in common may then be expected for any source location. 
Since galaxies are also found in underdense regions, the sim- 
ulations may in fact describe a plausible helium reionization 
scenario, since an appreciable fraction of galaxies may have 
undergone an AGN phase at some stage in their lives so that 
helium in some underdense regions may have been reionized 
by a local source residing within them. 

For brevity, the labels listed in Table [1] will be used 
henceforth to identify the various simulations with their 
combination of source location and source power-law turn- 
on redshift. 

All simulations, unless stated otherwise, were performed 
m a (25fe ~^Mpc)^ comoving volume. A ACDM model was 
assumed (jSpergel et al.l l2003). with parameters: h = 0.71, 
Qbh^ = 0.022, fi™. = 0.268 and 0,^ = 0.732, where ilt, fl^, 
and have their usual meanings as the contributions to 
Q from the gas, all matter, and the vacuum energy, respec- 
tively. We adopt a helium fraction by mass of y = 0.235. 

The initial density perturbations were created by dis- 
placing a uniform grid of 512"^ particles in the box using the 
Zel'dovich approximation. The initial power spectrum of the 
density fields was a COBiJ-normalised power law with index 
n = 0.97. The same initial conditions were used for all sim- 
ulations. Since there is no feedback from the RT to the PM 
code, all the runs have identical gas densities. The simula- 
tions were evolved to a redshift of 2. 

The reionization computations were carried out on a 
superimposed polar grid. Details of the scheme are provided 
in the Appendix. The finite travel time of light is included 
in the propagation of the radiation field from the source. Be- 
cause of the intensive computational demands required for 
accurate temperature determinations associated with reion- 
ization, the regions ionized are restricted to slabs 500h~^ kpc 
(comoving) thick and 25/2^''^h~^ Mpc (comoving) in radius. 
This produces a small amount of periodic overlap in the 
density field at the edges of the polar grid along its cardi- 
nal directions, necessary to maximise the usage of the cubic 
simulation volume. 
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Figure 3. The evolution of azimuthally-averaged line-of-sight gas temperatures around a central radiation source. The left column 
corresponds to a source placed at a density maximum while for the right column the source is placed at the centre of a void. From top 
to bottom, the source evolves from a starburst spectrum to QSO (power-law) over the redshift intervals 5.5 > z > 4.5, 4.5 > z > 3.5 and 
3.5 > 2 > 2.5, respectively. The different curves correspond to constant comoving impact parameters for the lines of sight of f) = 0.05 (solid 
line), 0.15 (long-dashed line) and 0.25 (short— dashed line), in units of the polar grid diameter. The temperatures show a monotonically 
increasing trend with impact parameter after the Hem I-fronts pass. The polar grid radius is 25/2^/^ Mpc (comoving). 
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Figure 4. The evolution of the azimuthally-averaged hne-of-sight median gas temperature for different over-density cuts, for run HD4.5. 
The different curves correspond to projected separations of b = 0.05 (solid line), 0.15 (long-dashed line) and 0.25 (short-dashed line), in 
units of the polar grid diameter. 
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Figure 5. The evolution of tlie relation between gas density and 
temperature in the vicinity of a source placed at a density peak. 
The dependence is split between the inner half radius of the po- 
lar grid (top panel) and outer half radius (bottom panel). The 
contour levels correspond to probability density contours per log- 
arithmic intervals in temperature and overdensity of 10""'^, lO"'^ 
and IQl-^. 



Figure 6. The evolution of the relation between gas density and 
temperature in the vicinity of a source placed at a density mini- 
mum. The dependence is split between the inner half radius of the 
polar grid (top panel) and outer half radius (bottom panel). The 
contour levels correspond to probability density contours per log- 
arithmic intervals in temperature and overdensity of 10""'^, lO"'^ 
and lO^-^. 



3 SIMULATION RESULTS 

3.1 Physical properties of the IGM 

We focus on the thermal state of the gas surrounding the 
central source before and after the Hem I- front passes. The 
evolution of the temperature in a slice through the simula- 
tion is shown in Figure [T] The role played by dense struc- 
tures surrounding the source in restraining the growth of the 
Hii and Hem I- fronts is apparent in the highly anisotropic 
temperature patterns. For comparison, the density field at 
z — 2.6 is shown in Figure[2l (Since the overdensities change 
very little over the redshift interval 2.6 < z < 4.6, only the 
image at z = 2.6 is shown.) The gas temperature increases 
sharply behind the Hem I-fronts, then relaxes to lower val- 
ues as the gas achieves equilibrium between photoionization 
heating and atomic radiative losses. As a consequence, the 
highest temperatures tend to lie furthest from the source in 
regions behind the Hem I-front. Thus, whilst the gas tem- 
perature is increased within the Hem region surrounding the 
QSO compared with the gas temperature outside, the tem- 
perature of the gas within the Hem region decreases towards 
the QSO itself. As will be discussed in ^below, another im- 
portant contributing factor to the trend is the hardening of 



the radiation field by the intervening gas between the source 
and the more distant regions. 

Figure |3] shows azimuthal averages of the median gas 
temperature along lines of sight drawn around the source 
with varying impact parameter. Once the power-law source 
turns on, a steep rise in the temperature with decreasing red- 
shift is found. The rise occurs over a much narrower redshift 
interval than the transition redshift interval of the source 
spectrum as the source becomes sufficiently hard to pho- 
toionize Hen. The associated photoionization heating occurs 
over a redshift interval of approximately Az ~ 0.4, with the 
temperature peaking at the centre of the transition redshift 
interval. The median averages show a trend of increasing 
temperature with increasing impact parameter of the lines 
of sight after helium is fully ionized. The median temper- 
atures approach convergence with time as photoionization 
heating comes into equilibrium with atomic radiation losses 
and adiabatic expansion cooling. 

Similar trends are found for both underdense and over- 
dense regions, as shown in Figure|4l with peak temperatures 
of r f» 20 X 10'^ K after Hen reionization, diminishing to 
r fa 10" K by 2 = 2. (Since the projected radh b > 0.25 
begin to sample the rising density field in the next periodic 
image of the simulation volume, subsequent plots are trun- 
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cated to 6 < 0.25.) At high overdensities, however, the trend 
is less clear and appears to reverse after helium reionization 
as a function of impact parameter, as shown by the panel 
for p/ (p) > 3. This density threshold corresponds to systems 
with Hi Lya optical depths of about unity, approximately 
the minimum for which meta l abso rption systems have been 
detected (eg, lAguirre et al.l (|2008h ). For this density cut, 
after the initial rise the temperature settles to roughly con- 
stant values depending weakly on the impact parameter of 
the line of sight, but now with lower temperatures for the 
larger impact parameters. The median temperatures are also 
found to be considerably higher than for most of the gas. 
The peak median temperature following full helium reion- 
ization reaches values of T « 30 — 40 x 10"^ K, decreasing to 
20 — 30 X 10^ K as thermal equilibrium between photoion- 
ization heating and atomic radiative losses is established. 

The trend of temperature with impact parameter is 
thus not a consequence of a narrow equation of state re- 
lating the gas temperature to the gas density. The post- 
ionization gas in fact shows a wide range of temperatures 
at a given overdensit y in reionization computations includ- 
ing ra,diative transfer llBolton et al.ll200i : iTittlev fc MeiksinI 
l2007l : lBohon et al.ll2009l )~ The dependence between gas tem- 
perature and density is shown in Figure [S] The relation 
has been split into two regions: an inner region within a 
half radius of the polar grid (top panel), and the comple- 
mentary outer region (bottom panel). At z = 4.4, just af- 
ter the power-law source turns on, the temperatures begin 
to elevate in the inner region as Hen is photoionized. The 
gas temperature and density obeys an approximate poly- 
tropic relation. This arises in underdense regions from a 
balance between photoionization heating and adiabatic ex- 
pansion cooling l|Hui fc GnedinI Il997l ). At densities above 
riH ~ 10~* cm~^ (corresponding to p/{p) > 3 — 10 for 
3 < 2: < 4.5), however, the gas density is sufflciently high 
that there is time for thermal balance to be established be- 
tween photoionization heating and atomic radiative cooling, 
which results in a trend of dec reasing temperature with in- 
creasing density (lMei ksir]|l994l '). This turnover is present in 
the outer region (lower panel). The extension of the poly- 
tropic relation to higher overdensities in the inner region 
reflects the role of adiabatic compression heating within col- 
lapsing structures, with T ~ p^^^ and reaching peak tem- 
peratures of mainly 3 — 5 x 10* K, as expected for colla psing 
struc tures on scales of 100 kpc (proper) and smaller (MeiksinI 
1 19941 ). The temperature-density slope is slightly shallower 
than the adiabatic relation (Figure [5|, as radiative cooling 
limits the temperatures in overdense structures. At virial 
densities (p/(p) ^ 200), for which the gas comes to rest, 
atomic processes lower the temperature, creating a down- 
ward trend in the relationship between temperature and 
density (upper panel of Figure [5]) . The highest temperature 
regions (T > 40000 K) in Figure [T] align with the densest 
structures in Figure [2J with the highest temperatures in the 
dense structure at the centre of the grid. 

By z = 4.0, the Hem I-front has reached the edge of the 
simulation volume. The temperature-density relation has re- 
laxed to a new polytropic relation in the inner region as adi- 
abatic expansion cooling lowers the temperature of the un- 
derdense gas with time (upper panel). The turnover at high 
densities is consistent with the dominance of atomic pro- 
cesses in establishing the temperature-density relation, as 



discussed above. By comparison, in the outer region, where 
the Hen was more recently photoionized, the temperature- 
density relation shows a pronounced flattening, gradually 
settling to a polytropic relation with a slightly steeper gra- 
dient with decreasing redshift as adiabatic cooling lowers 
the temperature of the underdense gas. In both regions, a 
wide spread in the relationship is found. The high median 
temperatures (T > 40000 K), in the inner region at high 
densities are not reached in the outer regions because of the 
absence of massive halos comparable to that situated at the 
centre of the slice. 

Similar evolution is found for a source turning on in an 
underdense environment, as shown in Figure [S] In this case, 
the inner region heats up more quickly after the source turns 
on, since the Hem ionization front is able to propagate more 
quickly. The faster expansion speed in the void surrounding 
the source results in an increased amount of adiabatic cool- 
ing, bringing the gas to lower temperatures in the vicinity 
of the source (top panel) . Since the gas does not reach high 
overdensities in the inner region, the high temperatures as- 
sociated with collapsing halos found for the source placed at 
a density peak do not occur. Further from the source (bot- 
tom panel), the temperature evolves similarly to the case 
with the source placed in a dense environment. 



3.2 Absorption line properties of the IGM 

Whilst isolated features in the H en Lyg forest ha ve been 
detected up to redshifts z < 3 (|Kriss et al.ll200lh , it be- 
comes increasingly difficult to detect individual Hen Lya 
absorption features at z > 3, as observations from space are 
required and the sources tend to be too faint at these wave- 
lengths to readily measure the spectra beca use of the hig h 
Hen Lya optical depths at these redshifts (|Meiksinll2007l ). 
We explore the possibility of detecting the effects of Hen 
reionization on the Hi Lya forest instead. The hydrogen 
ionization fractions were recomputed allowing for a back- 
ground metagalactic hydrogen ionization rate in addition to 
that of the QSO of, in units of lO"^'^ sT^ per neutral hydro- 
gen atom, 



0.8 



HI, -12 = 



4 
1 + z 



0.386 



6.5 

1 + z 



2 < z < 5.5 
5.5 < z < 6, 



(1) 



in go od a greement wit h re cent e sti mates (iMeiksin fc Whit3 



'2004'; 'Tvt ler et al. 
Fauchor- Gigu ere et al. 



l2004l : 
l2008l ). 



iBolton et all |2005| : 



Synthetic spectra were 
generated by casting lines of sight to fictitious background 
sources azimuthally distributed about the central ionization 
source. The spectra were computed from the modified Hi 
fractions using Eq. [T] and the gas density, temperature 
and peculiar velocity from the s imula tions, following the 
procedure in I Meiksin fc Whit"3 (|200ll ). Spectral regions 
were examined within ± 1000 km s"'^ of the closest point 
to the QSO along a neigbouring line of sight, a range 
sufficiently broad to produce a large number of absorption 
features whilst still within the region of influence of the 
central QSO. No corrections for the finite travel time of 
light were made in constructing the spectra except for the 
time delay between the source at any given position within 
the grid, which is small. Corrections are important for an 
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ionization front travelling at close to the speed of light, as 
it will near the source when the source first turns on. In 
principle, any comparison with observations could correct 
for the resulting time-delay effects in the observed spectra. 
The effects, however, are small over the scales presented. 
By Az = 0.1 less than the turn-on redshift of the QSO, the 
Hem I- front is expanding at less than 0.1c. The QSO ionizes 
most of the polar grid within Az = 0.3 of the turn-on 
redshift, with only highly subluminally expanding patches 
of Hem remaining in dense regions or regions obscured 
from the source by dense intervening lumps of gas. The 
spectra presented are generally drawn from the simulation 
either before or after the Hem I-front passes, except for a 
brief interval of width Az ~ 0.3 during which the I-front is 
moving sufficiently subluminally that finite speed of light 
corrections are at the 10 per cent, level or smaller, and so 
are not included. 

A Voigt-profile absorption line analysis was performed 
on the spectra using AUTOVP jPave et al.1 liggTl ). modified 
with additional error controls to render the code more ro- 
bust. Azimuthal averages of the median values of the result- 
ing Doppler parameters are shown in Figure [T] The Doppler 
parameters increase as the hard photoionization source turns 
on, reflecting the behaviour of the median gas temperature 
in Figure [3] although with a somewhat gentler rise. 

There is, however, a remarkable difference in the trend 
with projected separation of the line of sight. Prior to the on- 
set of the hard photoionization spectrum and the subsequent 
full reionization of helium, the gas temperature shows little 
variation with impact parameter from the central source. 
The small trend indicated in the figure is consistent with the 
trend with density as a function of offset, as shown in Fig- 
ure [S] After helium becomes fully ionized, the temperature 
increases systematically with distance from the source. By 
contrast, the Doppler parameter shows a strong trend of de- 
creasing value with increasing projected separation from the 
central source prior to the onset of the hard photoionization 
spectrum. The trend, along with a gradual decrease with de- 
creasing redshift, is apparent in the high density runs, par- 
ticularly HD3.5 and HD4.5, as well as in all the low density 
runs, as shown in Figure [T] 

After helium is fully reionized, the trends with impact 
parameter and redshift gradually disappear. In fact, regard- 
less of the onset of helium reionization in the models, by 
z — 3 the median Doppler parameters all take on nearly the 
same value of about 25 km s~^ in t he high density run s, in 
good agreement with observations ijMeiksin et al.ll200ll ). In 
the low density runs, a slightly lower median value of about 
23 km is found. 

The near constancy of the median Doppler parameter 
with both redshift and impact parameter following helium 
reionization is a surprising result. Since the Doppler param- 
eter increases with gas temperature, trends similar to those 
found for the temperature may have been expected. The 
Doppler parameter, however, reflects the underlying pecu- 
liar motions of the gas as well as the temperature, although 
the relation is a complex one. The contribution of pecu- 
liar motions to the broadening of the absorption lines is 
not straightforward to isolate since the motions not only 
broaden the lines but displace them as well. An attempt 
to construct spectra without including the peculiar velocity 
field resulted in qualitatively different absorption structures 



with relative displacements of line centres and substantially 
different blending of lines. Since absorption line fitting is 
a non-linear process, very different absorption lines result 
with widths that reflect differences in the deblending of the 
absorption features in addition to the absence of peculiar 
motions. The total amount of absorption changes as well, 
making it unclear how to renormalise the spectra for a fair 
comparison with the correctly computed spectra including 
the peculiar velocities. 

Some insight into the role of peculiar motions in broad- 
ening the absorption lines, however, may be made by exam- 
ining the structure of the peculiar velocity field. The line-of- 
sight peculiar velocity difference across a 100 kpc (proper) 
separation, correspon ding to the typical t hickness of the ab- 
sorption systems (eg, IZhang et al.l (|l998l )). is found to de- 
crease with decreasing redshift, as well as to decrease the 
more distant from the source, as shown for case HD4.5 in 
Figure [51 Prior to helium reionization, when the gas temper- 
atures are relatively low, the Doppler parameters are com- 
parable to the peculiar motion differences and follow the 
same trend of decreasing value with increasing impact pa- 
rameter. Once helium is fully ionized, the Doppler param- 
eters become nearly constant both with redshift and with 
impact parameter. This suggests they are no longer probing 
the mean evolution in temperature and peculiar motion of 
the gas, as weighted by volume, but rather the temperature 
in denser regions which shows little dependence on either 
redshift or distance from the source, as shown in Figure |4] 
for regions with p/{p) > 3. It is these overdense structures 
which correspond to systems with Lya optical depths near 
unity and above, which dominate the absorption in the spec- 
tra at z < 4, in contrast to the lower overdensity structures 
that dominate the absorption a t higher redshifts b ecause of 
their higher absolute densities l|Zhang et al.lll993 ). 

Similar trends are found for a source placed in the cen- 
tre of an underdense region. Prior to full helium reioniza- 
tion, the Doppler parameter for absorption features along 
neighbouring lines of sight increases with decreasing impact 
parameter of the line of sight (Figure [7|. This is found to 
correlate well with an increased difference in the line of sight 
peculiar velocity over the scale of 100 kpc for smaller impact 
parameters. At such small impact parameters, nearly the full 
peculiar expansion velocity of the underdense region will be 
probed, while at larger impact parameters the contribution 
of the expansion velocity to the line broadening will be re- 
duced by the decrease in the projection along the line of 
sight. After reionization, the Doppler parameters are again 
nearly independent of redshift and impact parameter, sug- 
gesting again that they are probing overdense structures for 
which the gas temperature has little dependence on redshift 
or distance from the source. 

The evolution in the Hi column density and Doppler 
parameter joint distribution is shown in Figure [9] for model 
HD4.5 over the redshift range 4.6 > z > 4.0. The onset 
of Hen reionization has two effects. It increases the median 
Doppler parameter, while narrowing the width of the Hi 
column density distribution as the recombination rate is re- 
duced due to the increase in the gas temperature. This en- 
hances the reduction in the number density of high column 
density systems along a line of sight due to the decrease 
in the absolut e density of the ga s as a result of cosmologi- 
cal expansion JZhang et al.lll99g 'l. A clear envelope of mini- 
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Figure 7. The evolution of the azimuthally-averaged Doppler parameters from absorption line analyses of synthesized spectra. The 
different curves correspond to projected separations of b = 0.05 (solid line), 0.15 (long-dashed line) and 0.25 (short— dashed line), in units 
of the polar grid diameter. 
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Figure 8. The evolution of the azimuthally-averaged line-of-sight median gas temperature, over-density, peculiar velocity difference (over 
100 kpc separations), and Doppler parameters from absorption line analyses to synthesized spectra, for run HD4.5. The different curves 
correspond to projected separations of 6 = 0.05 (solid line), 0.15 (long-dashed line) and 0.25 (short-dashed line), in units of the polar 
grid diameter. 
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Figure 9. The evolution of the Hi column density (in units of cm and Doppler parameter joint distribution for mock spectra drawn 
from model HD4.5. 

mum Doppler parameters, increasing with A^hi, is produced reionization, consistent with the increase in gas tempera- 
both before and after full helium reionization. For systems tures. 
with AThi > 10^*cm~^, the envelope is found to rise from 
&min « 7km to 6min w 10 — 15km s~^ following Hen 
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Figure 10. The evolution of the rate of heat input per Hell ionization surrounding a central radiation source which evolves from a 
starburst spectrum to QSO (power-law) over the redshift interval 4.5 > z > 3.5. The source is at a density maximum in the simulation 
volume. The evolution in the rates reflects the changing hardness of the incident ionising radiation field as filtered by the IGM. 




Figure 11. The evolution of the Hem ionization fraction surrounding a central radiation source which evolves from a starburst spectrum 
to QSO (power-law) over the redshift interval 4.5 > z > 3.5. The source is at a density maximum in the simulation volume. The images 
show the rapid increase in Hem ionization once the QSO turns on. A double greyscale is used to capture the dynamic range of the Hem 
ionization fraction, one ranging from 0.1 (black) to 1 (white) in the innermost regions near the source, and the other ranging from 10~^ 
(black) to 0.1 (white) in the outer regions in the panels a,t z = 4.4 and 4.0. 
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4 DISCUSSION 

Much of the temperature structure of the IGM after he- 
hum reionization may be accounted for by the hardening of 
the radiation field as filtered through the IGM during Hen 
reionization. Approximating the spectrum of the radiation 
field locally as a power law = JHeii(i^/i'Heii)~", where 
JheII is the angle-averaged intensity at the Hen photoelec- 
tric threshold frequency fHoii, the ionization and heating 
rates per Hen ion are, respectively, 

^Heii = ■^HciiCTHoii / dv — 



JHellCTHell 

/ip(3 + a)' 



(2) 



and 

CHeii 



JhcIICHc 



/•oo 



dv 



Z^HcII 



J^HcII 



JhcIio-HcII^'HcIi[(2 4- a) (3 -I- a) 



(3) 



where the photoelectric cross-section for ionising Hen is ap- 
proximated as (THcii(i^/i'Hcii)~^, and h-p is the Planck con- 
stant. This then gives for the expected energy injected into 
the gas per ionization 



^He II 'ipi^Hcii 



Hen 



2 + a 



4.36 X 10~^**(1-|-q/2)"^ J. 



(4) 



Deviations from this will occur due to spectral structure 
within the ionization front, especially as the mean free path 
of high energy photons is long compared with that of pho- 
tons just above the photoelectric threshold, resulting in a 
hardening of the spectrum through the I-front. The rela- 
tion serves to provide an approximate estimate for the ex- 
pected amount of energy injection and the sensitivity to the 
hardness of the effective spectral index within the ioniza- 
tion front, which may even achieve negative values (a < 0) 
in regions for which more high energy photons are able to 
penetrate than low. 

The evolution of the energy per ionization is shown in 
Figure [TOl At early times after the source just turns on (the 
panels for z = 4.4 and z — 4.0), heating rates per ionization 
more than an order of magnitude greater than Eq. Q for 
a = 0.5 (the source spectral index) are found. The helium 
in these regions is still largely in the form of Hei, with the 
Hem region just beginning to emerge from the source at 
z = 4.4, and not quite reaching the edge of the polar grid 
hy z — 4.00, as shown in Figure [TT] 

Whilst the Hen ionization fraction a;Hon evolves rapidly 
within the ionization front, it is found that the heating rate 
per ionization is nearly constant for a given value of xhcU- 
This is illustrated in Figure [121 The energy per ionization is 
well-described by 



10^ 



10^ 



-17.37+t/5+t^/2 



-16.85+t/10+(t-l)"72 



:t< 0.9475 



:t> 0.9475 



(5) 



where t = XHcH + Sa^Hci. The right region corresponds to gas 
that is ionising from Hei to Hen, whilst the left region shows 
the level of Hen as the helium becomes fully ionized. The 
bridge between the regions indicates the rapid depletion of 
Hei as it is converted into Hen. As the Hen is ionized into 




Figure 12. The energy per Hen ionization as a function of 
the Hei and Hen ionization fractions, xjjcl and xhcIIi respec- 
tively, forming narrowly-defined contours. The dependence is 
well-parametrised as a function of XHell + 2a;Hel (thin line). The 
contour levels are spaced by a decade. The left contours corre- 
spond to XHel = 0. The data are for z = 4.0. 



Hem, the energy per ionization rapidly declines as the effec- 
tive ionising spectrum softens towards the intrinsic spectral 
shape of the source. 

The expected boost in the gas temperature, allowing for 
the injected energy to be shared by all the species present, 
is 

2 1 riHell , 



AT = 



3fcB 2.23 nn 
16800 SHcii 



-hpi^cii 
10-" J 



K, 



(6) 



where fca is the Boltzmann constant. Comparison with Fig- 
ure [To] shows that the rises in the median temperature in 
Figure |3] following full helium ionization are of the magni- 
tudes expected, although the actual local increase may be 
much larger, depending on the evolution of the Hen ioniza- 
tion fraction a;Hoii within the ionization front. The gas will 
also rapidly cool in dense regions due to radiative losses, and 
in underdense regions due to adiabatic expansion cooling. By 
z = 4.00, the highest heating rates are confined to the pe- 
riphery of the polar grid (Figure [TO)) . where pockets of a;jjgjj 
continue to survive in shadows cast by overdense clumps 
blocking the source, as shown by a comparison between Fig- 
ures [11] and [2] In the inner regions, where the helium is fully 
ionized, the gas is able to cool through radiative losses and 
the temperature lowers, leaving the peripheral regions at rel- 
atively higher temperatures. This is the origin of the trend 
of increasing temperature with impact parameter shown in 
Figure [3] High heating rates are found to persist down to 
z = 3.60 even after most of the helium is fully ionized, and 
indeed increase by z — 3.20 as dense clumps continue to 
gather near the centrally placed source further filtering the 
radiation field, demonstrating that radiative transfer effects 
may continue to affect the temperature of the gas in distant 
regions even well after helium reionization has completed. 



5 CONCLUSIONS 

Using a coupled radiative-transfer and Particle-Mesh N- 
body code, we study the effects of helium reionization by 
a QSO source on the surrounding IGM. We consider QSO 
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source turn-on redshifts ranging from z — 3.5 to z — 5.5, 
with the QSO spectrum modelled as a power law. In the 
overdense environment expected for a QSO, the crossing of 
the Hem ionization front boosts the temperature behind the 
front to median values of T ~ 20 — 30 x 10'^ K, gradually 
relaxing to T « 10 x 10^ K by z = 2 as thermal equilibrium 
between photoionization heating and atomic radiative cool- 
ing processes and adiabatic expansion cooling is established. 
The peak temperature tends to be higher for earlier turn-on 
redshifts. A "thermal proximity effect" results with a height- 
ened IGM temperature within the Hem region produced. An 
identical source placed in an underdense environment shows 
similar behaviour. 

The median gas temperature surrounding the QSO, 
however, is not uniform; it shows an increasing trend towards 
the Hem ionization front, where spectral filtering by the in- 
tervening gas still hardens the radiation field from the source 
in patches and where helium was most recently ionized. We 
generate synthetic spectra to assess the possibility of detect- 
ing the trend in the absorption lines along the lines of sight 
to background neighbouring QSOs. The spectral signature 
of the post-helium ionization temperature trends nearly van- 
ishes in the Hi Doppler parameters, which show both little 
evolution and little dependence on distance from the source. 
For a source placed in an overdense region, a median Doppler 
parameter of 25 km s~^ at z = 3 is found, independent of 
impact parameter and the time of Hen reionization. A some- 
what lower median Doppler parameter of 23 km s~^ is found 
for a source placed in the centre of a void. The median values 
are in good agreement with the measured median Doppler 
parameter. This is in contrast to simulations without radia- 
tive transfer, for w hich the value is und erpredicted by a few 
to several km s~^ jMeiksin et al.|[200ll ). Late Hen reioniza- 
tion appears to account for the higher values. The absence of 
a dependence of the median Doppler parameter on redshift 
and impact parameter may be attributed to the absence of 
any temperature trend with redshift and impact parameter 
of the denser regions which dominate the Lya absorption 
spectra. By contrast, a trend is found prior to full helium 
reionization, with an increase in the Doppler parameters for 
a smaller impact parameter of the line of sight, reflecting 
similar behaviour in the peculiar velocity field. Thus the dis- 
appearance of a trend with impact parameter may indicate 
the onset of full helium reionization. 

A near polytropic relation is obtained between the gas 
temperature and density with a slope that is close to adia- 
batic, extending to overdensities of ~ 10^ in the vicinity of 
the source placed in a massive halo. The temperatures found, 
approaching 10^ K, are those expected for accretion onto ha- 
los and sheets, with radiative cooling lowering the temper- 
ature at higher overdensities. In regions without a massive 
halo, the turnover in the temperature occurs at overdensi- 
ties of only a few, again as expected from radiative cooling 
losses. No inverted temperature-density relation is found in 
the underdense gas. Whilst the temperature-density rela- 
tion flattens soon after Hen reionization, the temperature 
in the underdense gas follows a general trend of increasing 
with density, although with a wide spread in values. 

The detection of metal absorption systems would al- 
low a more direct probe of the gas temperature. Because 
the peculiar velocity contribution to the velocity width of 
an absorption line would be identical for features arising 



from different elements, while the thermal contribution (in 
quadrature) decreases like the mass of the ion, it is possible 
to separate these two contributions provided absorption fea- 
tures due to two or more elements from the same structure 
are detected. Such temperatures have been measured for a 
few systems, although the corresponding Hi column densi- 
ties tend to be high, typically exceeding 10^^ cm~ '^, corre- 
sponding to overdensities of a few to several (eg, iMeiksinI 
(2007)). The temperature trend as a function of impact pa- 
rameter is less clearly defined for such overdensities, with a 
weak trend of higher values for smaller impact parameters. 
None the less, a spread with impact parameter could indi- 
cate the onset of full helium reionization, particularly if very 
high temperatures exceeding 30 x 10"^ K are found at z < 4. 

It would be very interesting to use metal absorption 
features arising in the very diffuse IGM, corresponding to 
Hi column densities of 10^'^ cm^^ and lower, to probe the 
temperature structure of the IGM, both before and after 
full helium reionization. There are currently no detections 
of individual metal absorption features in such low density 
diffuse systems, but neither do observations preclude them. 
Exploiting them for temperature measurements, however, 
may need to await spectroscopy on a much larger future 
generation of telescopes such as the Thirty Meter Telescope 
or an Extremely Large Telescope. 
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APPENDIX A: RADIATIVE TRANSFER ON A 
POLAR GRID 

The equation of radiative transfer is solved along radii 
emanating from the source using a photon-conserving 
probablistic method jAbel et al.l Il999l : iBolton et all l2004l : 
iTittlev fc MeiksinI lioOTi r The radiative transfer was per- 
formed on a polar grid with the start of each line of sight 
placed at the source. The polar grid is divided into sectors 
bounded by intervals in radius and angle. Since the optical 
depth must be computed from the source to each polar cell, 
a density field must be defined on the polar grid in order 
to compute the column densities along the radii from the 
source. We considered three methods to generate a smooth 
density field matching that of the Cartesian grid. We de- 
scribe the final method adopted first, but briefly describe 
two alternative methods as well that were found inadequate 
for information only should any of these methods be at- 
tempted. 

In the method used, the PM mass distribution is ap- 
portioned to the cells in the polar grid by dividing each par- 
ticle into 500 sub-particles, spread over a comoving radius 
of 0.25 h~^ Mpc, and assigning the contributions to the po- 
lar grid cells in which each sub-particle lies. The sub-particle 
density distribution is cone-shaped, declining like 1/r in 2D, 
so that the number per ring of fixed width dr is constant. 
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(The corresponding distribution in 3D would decline like 
for an equal number per shell.) The mass density on 
the polar grid is then computed as the total mass in each 
cell divided by the cell volume, which increeises with radius. 
The method was found to produce a smooth rendition of the 
density field without discontinuities and in good agreement 
with the original density field measured on a Cartesian grid. 

Wc considered two other methods that, though seem- 
ingly straightforward, produced undesirable artefacts. In the 
first, the density field defined on the Cartesian grid was re- 
mapped onto the polar grid by computing the total mass 
contributed by the Cartesian grid cells that overlapped a 
given polar cell. Because of the irregular shapes of the over- 
lap regions, it was found impractical to compute the regions 
of overlap directly. We instead computed the overlap regions 
using Monte Carlo integration by populating the Cartesian 
grid uniformly with a large number of randomly distributed 
points in order to compute the fraction of overlap with in- 
dividual polar cells. Each Cartesian grid cell was then as- 
signed a weight for each polar cell according to the amount 
of overlap. Once computed, the weights were stored, as they 
need not be computed again for a given source placement. 
The method, however, was found to produce discontinuities 
in the density field if the source is placed near a Carte- 
sian cell boundary. The resulting density field on the polar 
grid was also found to be excessively smoothed compared 
with the density field defined on a Cartesian grid. A second 
approach was to compute the density field directly on the 
polar grid using the N-body particles. A nearest grid point 
method based on counting the particles in individual polar 
cells fails, particularly near the centre of the grid, because 
the cells are so small that large Poisson fluctuations in the 
density field Hr(> produced, including empty cells with zero 
density. Enlargening the cells near the centre results in poor 
resolution and over-smoothing. 
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